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Abstract. In this work we show that 3-3-1 model with right-handed neutrinos has a natural 
qq ■ weakly interacting massive particle (WIMP) dark mater candidate. It is a complex scalar 

with mass of order of some hundreds of GeV which carries two units of lepton number, a 
scalar bilepton. This makes it a very peculiar WIMP, very distinct from Supersymmetric or 
C"| . Extra-dimension candidates. Besides, although we have to make some reasonable assumptions 

concerning the several parameters in the model, no fine tunning is required in order to get the 
O correct dark matter abundance. We also analyze the prospects for WIMP direct detection by 

considering recent and projected sensitivities for WIMP-nucleon elastic cross section from 
CDMS and XENON Collaborations. 
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1. Introduction 



CN 

d 

^ ; The problem of matter density in the Universe seems to be one of the most intriguing and 
exciting subjects in modern Physics. The growing refinement achieved in cosmological data 
leaves no doubt about a dark component in the observed mass density, constituting roughly 
22% of all energy density acording to the three year run of WMAP \M ■ This yet unknown 
component has to be non-baryonic, its interaction with the electroweak Standard Model (SM) 
particles should be negligible and it has to be cold, i.e., non-relativistic at the time it decouples 
from the radiation bath, the so called Cold Dark Matter (CDM). From the theoretical side, 
there are some proposals to explain the CDM in the context of Particle Physics models 
(see Ref. []2|, [3], |], |J, and references therein for a review of the subject). Among them 
there are models which present natural candidates to play this role, the weakly interacting 
massive particles (WIMP)'s, with mass ranging from approximately 1 GeV to 1 TeV. These 
WIMP's are nice candidates because their masses are in the GeV realm, turning them cold at 
decoupling, and mainly because their weakly interacting aspect not only yields a thermally 
averaged annihilation cross section of order of weak interactions, leading to the expected order 
of magnitude to CDM abundance, but also coincides with the scale of Particle Physics models 
to be probed at the Large Hadron Collider (LHC) that finds itself at the final stage to start its 
running phase [|7|]. It also presents the possibility of being seen in direct detection experiments 
since its massiveness would imply an observable recoil of nuclei in elastic collisions [EL Rp. 
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The most promising scenarios where such WIMP's can be present in the particle 
spectrum are Supersymmetry (SUSY) and Extra Dimensions models [Q, ^ |5p. All these 
models dispose of some kind of discrete symmetry in order to stabilize their CDM candidates. 
Also, they have to be realized at the electroweak scale so that their new particles are potentially 
good candidates for CDM. Although such models may represent the greatest expectations 
in Particle Physics for the Physics at TeV scale to be probed by LHC, the absence of any 
experimental evidence to support these models allows us to work with alternatives. One 
of such alternative routes concerns the enlargement of the gauge symmetry group from 
SUc{3) ® SUl(2) (g) Uy(l) to larger groups. In particular, there exists a simple extension 
of the SM gauge group to SU C (3) ® SU L (3) ® U x (l), the so called 3-3-1 model 0. 

This class of models is interesting for several reasons, not only because they mean a 
different scenario, but because they possess several nice features. For example, (i) the family 
problem is absent in this model since it demands that there be only three families of fermions 
when anomalies are canceled and asymptotic freedom is considered [|9|]; (ii) electric charge 
quantization is automatic [|T0[]; (iii) right-handed neutrinos can be part of the spectrum in 
some versions of the model [ |TT| , D"2| , |T3"1 ] and their tiny observed mass difference can be easily 
accommodated []14{]; (iv) axions and majorons are a natural outcome in some versions [|I3[], 
providing light particles which could also contribute to the problem of dark matter origin. 
Besides, it is possible that a custodial symmetry exists in these models which would make 
them indistinguishable from SM at low energy scales [p^||, and it would then be a strong rival 
to SM itself. 

In this work we will concentrate on the 3-3-1 version of the model with right handed 
neutrinos in the spectrum (3-3-lRHv) model [|T2"|]. The reason behind this choice relies on 
the fact that neutrinos mass is already a mandatory property that needs to be included in all 
reasonable extensions of SM. Besides, the model can be implemented with just three scalar 
triplets instead of including a sextet as in other versions, being economical in its content. 
However, the most appealing motivation to deal with 3-3-1 RHv to explain the origin of 
CDM is due to the possibility of having a candidate with a very distinct signature. Among 
its properties the model can be made lepton number conserving if some of its fields carry 
two units of lepton number which will be called bileptons. This peculiar property has many 
phenomenological implications. Namely, rare lepton decays can emerge, neutrinoless double 
beta decay is allowed, right handed neutrinos are going to appear as byproducts of heavy 
vector bileptons decays and so on. It is then automatic to ask if some of these additional 
bilepton fields can be a CDM candidate, once it is provided with a very specific quantum 
number appropriate to forbid its interaction with many of the electroweak fields. This would 
play a similar role as that played by the discrete symmetries in the competing models cited 
above. Moreover, as the sought candidate is merged in the exotic new effects just mentioned, 
their appearance in the coming collider experiments would represent an unquestionable 
evidence of our CDM candidate. 

What we are going to investigate in this work is the possible realization of this scenario 
in the 3-3-lRHv model for one of the bilepton scalars. We are interested in identifying this 
field and characterize it as a WIMP. This can be realized if it can be shown that it is stable 
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in the range of parameters for which the abundance is in accordance with the recent data of 
WMAP [p. It is important to stress that the natural perturbative scale for the 3-3-1 models 
is on the TeV scale [|T7|], which is suitable for obtaining a WIMP. In view of all this, it seems 
that such a possibility is as welcome as any previous attempt to explain the CDM content 
through SUSY or Extra dimensions, offering a completely new and distinct particle to do this 
job. It should be mentioned that other works exist in the literature trying to explain CDM in 



the context of 3-3-1 models Q18I], but their aim was to obtain a self-interacting dark matter to 
avoid excessively dense cores in the center of galaxies and clusters as well as excessive large 



number of halos within the local group when contrasted to observations Q19Q, which demands 
a light dark matter candidate. We do not pursue this approach here. On the contrary, we 
want to show that the 3-3-1 RHu model possesses a bilepton scalar that can play the role of 
a WIMP, which is the preferred candidate for CDM. Besides, we get this without the need to 
fine tune its couplings to small and unnatural values. 

This work is organized as follows. In Sec. [Z] we present the model with its content and 
interactions. In Sec. [3] we diagonalize the mass matrices for the particle spectrum, allowing 
us to characterize the new extra particles. In Sec. ^ we identify our WIMP with one of the 
neutral scalar bileptons of the model, checking its viability as CDM candidate by computing 
its abundance for a range of values of the parameters, which turns out to be natural in the sense 
we do not need to make any fine adjustment on the parameters. Then we study the prospects 
of direct detection for this WIMP. We conclude with Sec. ^. 



2. The model 



In the 3-3-lRHz/ model the leptons come in triplet and in singlet representations, 



aL 



e a L 



1,3,-1/3), e aR - (1,1,-1), 



(1) 



while in the quark sector, one generation comes in the triplet and the other two compose an 
anti-triplet representation with the following content, 
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(3,3,0),«iR -(3,1,2/3), 



(3,1,-1/3), d' m -(3,1,-1/3), 
fu 3 \ 



d 3 

V M 3 J 



(3,3, 1/3), u 3R -(3,1,2/3), 



(3,1,-1/3), u' 3R -(3,1,2/3) (2) 

where a = 1, 2, 3 refers to the three generations and the index i = 1, 2 is restricted to only 
two generations. The primed quarks are new heavy quarks with the usual electric charges. 
Actually, as we will see below, they are leptoquarks since besides baryon number they also 
carry two units of lepton number. 
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In order to generate the right masses for the known particles the model requires only 
three scalar triplets, namely, 



(3) 



with 7] and x Dom transforming as (1 , 3 , —1/3) and p transforming as (1 , 3 , 2/3). 

We assume the following discrete symmetry transformation for the fields in order to have 
a minimal model, 

(x , P , e aR , u a R , u' 3R , d' iR , Q 3L ) -> - (% , p , e aR , u aR , u' 3R , d' iR , Q 3L ) (4) 

This symmetry helps in avoiding the undesirable Dirac mass terms for the neutrinos and 
allows a minimal charge conjugation and parity (CP) conserving potential []20[], 

V(V, Pi X) = P^X 2 + ^V 2 + P 2 P P 2 + -^iX 4 + A 2 ?7 4 + ^3P 4 + 
A 4 (x t x)(^ t ^) + A 5 (x t x)(p t p) + K{rfv)(p ] p) + 

\ 7 (x^v)(v^x) + A8(x f p)(p f x) + Mv i p)(p t v) 

f 



V2 



e ijk ViPjXk 



H.c. 



(5) 



Also, the Yukawa sector can be written as, 

- C Y = fijQiLX*d'j R + hzQzLXu' 3R + g ia QiLV*daR 

+ h 3a Q 3L rju a R + g3aQ3Lpd a R + h ia Q iL p*u aR + G aa f aL pe aR + H.c. (6) 

With this set of Yukawa interactions all fermions, except the neutrinos, gain mass. In this 
model neutrinos masses are generated through effective dimension-five operators as shown in 
Ref. 0]. 

In the gauge sector, the model recovers the usual SM gauge bosons, , Z° , 7, and 
contains five additional vector bosons called V ± , U°, U ^ and Z' [12], with masses around 
hundreds of GeV. 

Considering the required properties for a dark matter candidate, it will be important to 
have in mind that some of the new particles carry two units of lepton number L and are known 
as bileptons, namely, 



Uv + , t/ t0 , u 3 , v 

This assignment is such that the lagrangian is lepton number conserving. Observe that this 
quantum nunber association limits the range of interactions available to the bileptons since the 
only gauge bosons that carry lepton number are the new ones, and their masses are at the TeV 
scale. This property will show itself useful when considering the WIMP stability in Sec. |]. 



J ° P + ) 



-2, L(y-,c/ Q ,^ 



X 5 X 



+2. 



(7) 



3. The mass eigenstates 

From the scalar content of the model Eq. (Q) we have five neutral scalars at our disposal. Two 
of them carry two units of lepton number, r]'° , x° an d if lepton number is conserved (as we 
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assume here) they do not develop vacuum expectation value (VEV). As for the remaining 
ones, since they do not carry lepton number they are free to develop nontrivial VEV's, 

V°, P°, X° -> 4|(^, AX ' + R ri,P,x' + «W)' ( g ) 

On substituting this expansion in the above potential, we obtain the following set of 
constraints, 

2 , \ 2 , ^4 2 , ^5 2 f v v V P n 
+ A i*V + + ~2 Vp ~ 2 "77 = ' 

/^^ + ^, + ^-|^ = 0. (9) 

It is reasonable to assume that / is of the order of v x >, the scale associated to 3-3-1 breaking 
ocurring at few TeV, while v v and v p are related to electroweak symmetry breaking. In what 
follows we assume, for simplicity, that the VEV's related to the ordinary vector bosons mass 
are the same, and use the convenient notation, v v = v p = v, and also / = v x >/2 = V/2. 

By construction, lepton number is conserved by the interactions of this model, implying 
that the neutral scalars (x° , do not mix with the other three x'° > ?7° , p°. Thus, in the 
basis (x° , we have the following mass matrix HJ, 

T (Ar + l/2)^ 2 Jj. (10) 

After we diagonalize this matrix, we obtain a zero mass scalar which is given by 

G= ^—x" + , 1 rf>, (11) 

recognized as the Goldstone boson eaten by the gauge bosons U° and {7°t. The other scalar, 

' + (12) 



is a heavy scalar with mass given by M| = |(A 7 + |)(f 2 + V 2 ). The assumption V >> v, 
allows us to say that the Goldstone bosons are mostly contained in the complex scalar G ~ x°> 
while the heavy scalar is mostly <p ~ ?? /0 - We are going to show in the next section that has 
all the appropriate features to be our WIMP candidate for dark matter. Meanwhile we perform 
the mass matrix diagonalization for the remaining neutral scalars as well as the charged ones. 
We do that because we need to have control of all masses in the model in order to guarantee 
the stability. 

Considering the expansion in Eq. @, we have the following mass matrix for the CP-even 
scalars in the basis (R x i , R v , R p ), 



I M -t- S y 2 { 2 8 ) y { 2 s)y. \ 
+ A 2 



V 



' A4 1\ D . 1 \ V 2 1 j_ \$ V 2 



2 8W A2 V 2 8 2 V 2 

{^5 _ l\v_ _I _|_ Ae v 2 1 I \ v 2 
\ \ 2 8 J V 8 + 2 V 2 8 + A3 V 2 I 



(13) 



% Here we are considering the CP even and CP odd scalars altogether, \° — Rx° + ^x" anc ' V' — + Hn' , 
since they possess the same mass matrices. 
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Its diagonalization leads to the following scalar eigenvectors, 

H ^ ^=(R V + R p ) , 
whose masses are respectively, 

v 2 

m 2 h o w Ai\/ 2 + — , 



1 



(4A 2 + 4X 3 )v 2 + V 2 l + Wl-8A f 



m 



H 



(4A 2 + 4X 3 )v 2 + V 2 1 - \ l - 8A e 



^ 2 



(14) 



Among these, the first and second scalars in the above equation are heavy while the third one 
is the lightest scalar of the theory, which we recognize as the standard Higgs boson. 

Regarding the CP-odd scalars, we have the following mass matrix in the basis 



V 2 



( 4 

V 2 



V 



V 

v 
V 



V ' 



1 / 



(15) 



This matrix has the following eigenvectors, 



i 2 ~ 



7° 



V2 



(16) 



where P{ and 7 2 correspond to null eigenvalues and 7g is the massive pseudo-scalar with 



! v 2 



(17) 



The pseudo-scalars, 7° and 7°, are identified as the Goldstones eaten by the neutral gauge 
bosons, the standard Z° and the Z' characteristic of the extended 3-3-1 gauge symmetry. 

Finally, let us consider the charged scalar mass matrices. Remember that two of the 
charged scalars carry two units of lepton number (charged scalar bileptons), namely, x~ , p'~, 
while the other two, if , p~, have zero lepton number. This means that the first two charged 
scalars do not mix with the last two. Thus we are going to have two 2x2 mass matrices. The 
first one, in the basis (x~ > p' ) takes the form, 
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The other mass matrix of charged scalars, in the basis (77 , p ), is given by 

Their diagonalization shows that two of the four eigenvalues are null and correspond to the 
eigenvectors, 

fcJ«^(»T-/0. (20) 
while the massive states are, 

h 3 » p _/ , 

Again, the massless charged scalars are Goldstones, h~[ is eaten by the charged gauge boson 
V~ (recall that both are bileptons), and is the Goldstone eaten by the standard gauge 
boson W~. The charged scalars remaining in the spectrum, (bilepton) and hj, possess the 
following masses, 

<-4(A 8 + i)(F 2 +, 2 ), 

2 V 2 o 

m h 4 - = -y + A9W • (22) 

This completes our scalar bosons analysis, assuring us of the correct number of Goldstones, 
the identification of the Higgs boson and the extra heavy scalars. We observe that we can 
choose the scalar coupling constants of order 0.10 to avoid any fine tunning. Also, this 
choice turns the potential stable and all eigenvalues of the mass matrices positive. Finally, 
to guarantee that <p be the lightest new particle in the spectrum we have only to assume that 
A 7 is negative, which does not alter the stability of the potential and proves to be a choice as 
good as a positive coupling. 

Bellow we present the vector bosons masses of 3-3-\RH v model, so that we can compare 
them with the scalar bosons. We assume the experimental values for the SM gauge boson Z, 
M z = 91.118 GeV, while 

2 1 2 2 
m w± = 2 9 v ' 

2 2 1 2/ta2 1 2\ 

m v± = m u0 = j9{V + v ) , 

m% =^r9 2 (V 2 + v 2 ), (23) 

where we have used the approximation V ^> v and neglected terms of order higher than 
0(v 2 /V 2 ). The above degeneracy for the vector bileptons is due to our simplifying choice 



v v = v p = V. 



Next we expose our reasons to choose (f> = r/ u as our candidate for the 3-3-1 WIMP. 

§ Some of these couplings can be adjusted around this value in order to yield a Higgs mass of 1 14 GeV or a 
little higher, but no fine tunning is required for that. 
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4. The 3-3-1 WIMP 

A good CDM candidate should be an electrically neutral particle, be stable and generate 
an energy density about 22% of the critical energy density in the universe. The 3-3-lRHz/ 
model particle which fulfils these criteria, as we next show, is the complex scalar 0, which is 
electrically neutral, can be stable if it is the lightest new particle in the spectrum and, for some 
range of the parameters, can be shown to possess the right abundance for CDM. 

Concerning stability, notice that it carries two units of lepton number L and since the 
model is lepton number conserving, its decay can proceed only through final states with the 
same total L. Then triple interactions can be generally cast as r tr ipi e = (0) x (A) x (£?), 
where A is any particle carrying two units of lepton number, as given in Eq. ([7]), and B is 
any particle possessing null lepton number. By assuring that particles carrying L= ±2 which 
interact with are heavier than 0, its stability is then secured. 

By considering the Yukawa interactions Eq. (^|), we see that does not couple to leptons, 
but does couple to the new heavy quarks, which also carry two units of lepton number, u', 
d! or s', plus an ordinary quark. The bilepton quarks are expected to be heavier than since 
their masses are proportional to V, being of order of few TeV, thus forbidding to decay 
into fermions. As for the vector bosons as final states, once the decay has to involve V ± or 
U° to conserve lepton number, it is forbidden if is lighter than V ± or U°. This will be 
the case for almost the whole range of values of the free parameters which lead to the correct 
abundance and Higgs mass, as long as A 7 < in the potential Eq. (Q). Also, from the potential 
Eq. (Q) and the lepton number assignment Eq. (0), all triple scalar interactions with involve 
simultaneously hf and hf, in order to guarantee lepton number and charge conservation. 
However, as can be seen from Sec. [3|, all scalars except the Higgs are heavier than given 
that A 7 < and the other couplings are of order of one [j], thus cannot decay into a pair of 
scalars. These are all interactions we need to consider to be sure that is stable. 

Next, when computing the thermal averaged cross section, we establish the set of 
parameters that realizes this scenario. 

4.1. Relic abundance 

Once we have identified which particle can be our WIMP CDM candidate, it is imperative to 
obtain the correct observed CDM abundance VLcdm- We consider that the WIMP is in thermal 
equilibrium with radiation in the early epochs until its rate of reaction becomes smaller than 
the rate of expansion of the Universe, i.e. it decouples from thermal bath, the so called freeze- 
out. In the case of WIMPs this happens when they are non-relativistic. The Boltzmann 
equation dictates the evolution of the particle number density n with the expanding Universe, 

— + 3Hn = - < av re i > (n 2 - n 2 eq ), (24) 

where, H is the Hubble parameter or expansion rate of the Universe, which can be written 
as H 2 = 87rp/3Mp Z for a flat Universe, < av re [ > is the thermal averaged cross section 

|| See also in Fig. || the mass of the particles which interact with (j> as a function of the free parameters relevant 
for this work. 
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for WIMP annihilation times the relative velocity, n eq is the particle number density at 
equilibrium and M Pt w 1 x 10 19 GeV is the Planck mass and p is the energy density of 
the Universe. Since WIMP is non-relativistic at the time of decoupling (T << M^), its 
equilibrium number density is, 

n eq = 9<t> [-^) e T ' ( 25 ) 

with the WIMP number of degrees of freedom. Next we follow the standard procedure 
derived in Ref. [0, |2T]] to obtain a solution to the above equation and then determine the 
WIMP abundance. By defining x = M^/T and using the non-relativistic approximation for 
the squared center of mass energy, s = AM"^ + M^v 2 , we expand the cross section till the first 
power in v 2 , resulting in the thermally averaged cross section, 

6b 

< av rd >^a-\ , (26) 

x 

where a and b are the model dependent parameters. The relic WIMP abundance in the context 
of a flat cosmological constant dominated Universe (ACDM model), can be expressed as, 

„ , 2 1.04 x 10 9 x F 

VLsh 2 « — — — - ,. N , (27) 

where i F is i = M^/T computed at the temperature of decoupling of dark matter from 
equilibrium, the freeze-out temperature given by, 

^45 g+ M^M Pi (a+f; 



Xf = In 



c(c + 2)i 



(28) 



8 27T 3 y/g^xT 

which in general is close to ip ~ 20. Also, g* is the number of degrees of freedom of 
relativistic particles in thermal equilibrium with the WIMP at freeze-out, which in our case is 
about g* = 183/2 = 91.5 and g^ = 2. The constant c of order of one is obtained by matching 
the late-time and early-time solutions for the abundance and, for our purposes, it is enough to 
take it as c = 1/2 since it has only a small effect in the logarithmic dependence of xf- 

We can solve Eq. [28] iteratively and plug the result in Eq. |27| in order to compare the 
WIMP abundance with the latest results of WMAP [[p which, according to the ACDM model, 
imposes the following bounds to the dark matter abundance at 2a level, 

0.096 < n DM h 2 < 0.122. (29) 

The processes contributing to < av rd > are mainly the annihilation into gauge bosons, 
4>4> — > W + W~ and (jxp —> ZZ depicted in Fig. [II the annihilation into Higgs boson, 
00 — > HH shown in Fig. 0, which is comparable to the gauge bosons contributions, and 
the annihilations into quarks 00 — > qq presented in Fig. [3], whose role is better discussed 
bellow. 

Our first step to calculate the thermally averaged cross section parameters, a and b 
appearing in the above equations, is to set the range of free parameters of the model. In 
order to proceed with this we make some simplifying assumptions since there are too many 
free parameters to deal with. A natural choice is to take the couplings not far from 0(1). 
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Figure 1. WIMP annihilation into Standard Model gauge bosons, 

6*6 ^> ZZ. 
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Figure 2. WIMP annihilation into Standard Model Higgs, (f>*(f) — > HH. 
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Figure 3. WIMP annihilation into Standard Model quarks, <f>*(j) — > qq. 



In this sense, and considering the relevant interactions (see the appendix) involved in these 
annihilation channels, we leave the couplings, Ai and A2, free and fix 

A3 = A4 = A5 = Ag = Ag = Ai , 

A 6 = - A 7 = A 2 . (30) 

This convenient set of parameters allows us to vary the Higgs and <\> masses as we wish and 
the only care we have to take is not to make the bilepton gauge bosons, V ± and U°, lighter 
than (f). As we will see next, this is the case if V > 1.2 TeV and A 2 > 0.078, which is true for 
all parameter space which is compatible with WMAP results for VLdm- 

Concerning the annihilation into quarks, we will assume that the Yukawa couplings are 
dominant for the quarks in the same family and much smaller otherwise. These assumptions 
allow us to neglect any mixing among quarks in this study. In this case, according to Eq.@, 
the relevant couplings are given in terms of down, strange and top mass, g n = \f2m d /v v , 
922 = y/2m s /v v and /i 33 = \Z2mt/v v , respectively. The WIMP annihilation into quarks are 
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restrict then to the following reactions, 

00 — > dd , 00 — > SS , 00 — > it . (31) 

It should be stressed that although the annihilation into top quark is the only one where the 
s channel can be dominant since it proceeds by an intermediate Higgs, the same is not true 
for annihilation into down or strange quarks. This happens because the exchange channels 
are not suppressed by the lightness of quark mass as happens for the s channel. Indeed, for 
the Yukawa couplings assumed here, the amplitudes for t and u channels are some orders 
of magnitude higher than the amplitude for the s channel in the process 00 — > dd and 
00 — > ss. Nevertheless, when we compare the magnitude of these contributions with the 
s channel amplitude for 00 — > it, they are negligible. For this reason we only take into 
account those diagrams in Fig. |3| for top quark annihilation. 



Summing up all these contributions to the WIMP annihilation flfl, we use equations (|26|), 

to see if there is any compatible range of 



( ]27[ ) and ( p8| ) and impose the bound in Eq 
parameters which gives the correct Dark Matter abundance. Our results are shown in Fig. |]. 
In this figure we plotted the abundance against the coupling A 2 for fixed Ai, and the other 




0.05 



0.1 



0.15 



0.2 



Figure 4. The WIMP abundance as a function of the coupling A2, with Ai = 0.7 and other 
couplings constrained to these ones as exposed in the text. From the right to the left, the green 
curve is for V = 3 TeV, the red one is for V = 2 TeV, the blue is for V = 1.3 TeV and 
the black curve is for V = 1 TeV. The narrow shaded band encloses the allowed region for 
0.096 < n DM h 2 < 0.122 according to WMAP. 

couplings given in Eq. (J30|). For the range of values we took for the parameters, all scalars 
(except the Higgs) are always heavier than the WIMP, as shown in Fig. |5][[]. Since the Higgs 
mass is above M H = 114 GeV only for Ai > 0.7 (for V = 2 TeV), we used this value to 
get the abundance curves once the results are very weakly dependent on Ai. From Fig. |4] we 
observe that the larger the value of V, the larger is the maximum allowed value for A 2 . 

The expressions for the a and b parameters in Eq. (^6|) are too lengthy and we omit them in this work, though 
they can be easily obtained through the interactions listed in the appendix. 

+ Those plots were taken for V — 2 TeV, but they possess a similar behavior for other values of V. 
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The curve for V = 1 TeV presents only one valid region for 0.049 < A 2 < 0.063 leading 
to an acceptable dark matter abundance. However, for this value of V the Higgs mass is below 
its experimental lower bound, Mh > 114 GeV, for Ai = 0.7. If we change this coupling to 
larger values, for example Ai = 0.9, the acceptable region is still discarded for V = 1 TeV. 
This situation changes when we increase V to a minimum value V = 1.3 TeV, bringing the 
model to a comfortable position concerning the Higgs mass, even for X 1 = 0.7. However, 
when we look to the WIMP mass, we conclude that also the V = 1.3 TeV has problems 
for values of A2 < 0.078, since in this regime the vector bileptons would be lighter than our 
neutral scalar bilepton (see Fig. |5[), jeopardizing its stability. This leaves only a tiny window 
for A2 in this case, 0.078 < A 2 < 0.081. For V = 2 TeV and above, this problem is absent, 
though the possible range for A 2 is still short. We stick with values for V above 1.3 TeV and 
not much bigger than 3 TeV, in order to not reach the non-perturbative regime of 3-3-1 model, 
which would be attained at V ss 4 TeV JT7|] . 
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Figure 5. Mass spectrum for the additional scalars and vector bosons of 3-3-lRHi/ model for 
fixed V — 2 TeV. In (a) we show the neutral heavy pseudo-scalar as a function of A2, for 
Ai = 0.7. In (b) we plot the charged scalar masses valid for all values of A2 as a function of 
Ai. In (c) we show the WIMP mass against the vector bileptons masses as a function of A2, 
for Ai = 0.7. Finally, in (d) we plot the Higgs mass for two values of Ai = 0.7 and 0.9 as a 
function of A2 . 

It is appropriate to say that the above results were obtained for the bilepton quarks 
degenerated in mass, an assumption made to simplify our calculations, but departure of 
this should not modify our qualitative results. The only important point to be considered 
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is that bilepton quarks should have a mass larger than 1.5 TeV in order to generate sufficient 
abundance. This is a reasonable assumption since they receive mass from the largest scale in 
the model, v' = V, which is in the TeV range. Larger values for the masses of bilepton quarks 
would just push the values of A 2 to a bit larger values, without compromising our conclusions. 

We should also remark that we limited our analysis to a very conservative scenario, where 
several parameters were fixed ad hoc, since they are almost free of constraints. This does not 
mean that the above results put severe constraints on the couplings of 3-3-lRHz/ model if it 
has to contain a WIMP. We only meant to show that even with a restricted set of possibilities, 
the 3-3-lRHu model offers a good WIMP candidate for CDM with no need of adjusting its 
parameters to extremely unnatural values. Besides, the 3-3-1 RHv WIMP has a very peculiar 
signature once it carries two units of lepton number and could be easily distinguished from 
other models like SUSY or Extra Dimensions in collider experiments. 

From the above results, see Figs. (Q) and (|5|), our scalar bilepton 3-3-lRHz/ model WIMP 
has a preferred mass around M$ ps 600 GeV, which is about six times greater than that 
of neutralinos in the SUSY preferred scenario 0], and about the same magnitude as the 
lightest neutral vector boson Kaluza- Klein first mode, B^\'m 5D Universal Extra Dimensions 
(UED) [0]. The relative small WIMP mass in the case of SUSY can possibly be attributed to 
the fact that neutralinos are Majorana fermions and the main contribution to its annihilation 
occurs through P-wave into fermions. This leads to small cross sections requiring a smaller 
WIMP mass in order to annihilate more efficiently and give the correct CDM abundance. As 
for the in UED, its cross section can be shown to be roughly temperature independent [0], 
which means it efficiently annihilates into light fermions through S-wave allowing larger 
WIMP masses. Also, if the lightest WIMP of UED is a spinless photon in 6D, annihilation 
into fermions is helicity suppressed and presents a preferred mass half the size of the B^ 
case (see the last paper on Ref. [0]). 

Since <p in 3-3-lRHz/ model is a scalar, it is not S-wave suppressed for annihilation into 
gauge and Higgs boson, but as in the case of spinless photon in 6D UED, its annihilation into 
light fermions is also helicity suppressed. Indeed, for low values of V, the main contribution 
to the thermally averaged cross section comes from the annihilation into gauge bosons and, 
for some range of parameters, also into a pair of top quarks. Actually, as we have seen above, 
V has to be bigger than 1 .3 TeV, given the parameter space we are considering in this work, for 
which top pair annihilation is as important as gauge bosons to produce efficient depletion of 
leading to the right CDM abundance. This happens because P-wave is enhanced since the top 
quark is heavy and we choose a rather strong (p — T — t coupling. But also because the bilepton 
quark T exchanged in t and u channels is not extremely heavy, we took M T = 1.5 TeV, 
otherwise this contribution would be negligible and enough annihilation would occur for 
smaller <p mass. This is in agreement with our above results for VL C dm, which would demand 
increasing values for A 2 as M T increases. Thus, it is possible that such a similarity between 
the size of WIMP mass in 3-3-lRHz/ model and the vector boson in 5D UED is an artifact of 
the peculiarities of the chosen parameters in the former model. 
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4.2. WIMP direct detection 



After concluding that our proposed WIMP possesses some region of the parameter space 
which is in agreement with data on CDM abundance, we should at least check if some 
detection is possible for our candidate, or if there is any contradiction with the experimental 
exclusion limits put by the latest results on WIMP detection. The most probable signal of a 
WIMP is expected to appear in direct detection experiments, which consists of measuring the 
recoil energy of nuclei when these are elastically scattered by a WIMP [Q]. The data can be 
translated to a cross section normalized to nucleon and are usually presented as a limit on this 
cross section. We are going to analyze then the chance of direct detection of our scalar bilepton 
WIMP through elastic scattering with nuclei using the CDMS and XENON collaboration 
data which are the most stringent current results on WIMP-nuclei elastic collision. We 



will also use the projected sensitivities for WIMP detection in future experiments Q22|, |23Q. 

Our WIMP interacts with nucleons through its couplings with quarks by exchanging a 
Higgs boson or a bilepton quark (see Fig. The S matrix amplitude for these processes are 



H 



d.sA 



D,S,T 



\ 4> 




D,S,T 



Figure 6. Diagrams that contribute to the WIMP-nucleon elastic scattering cross section 
through WIMP-quark interactions. The first graph is the t channel Higgs exchange, while 
the remaining ones are s and t channel bilepton quarks exchange. 



given by, 



M a 



Mr 



Gf m i fO\ i \ \ — 



— I 



9, 



(Ms + E q ) 2 

n 2 



M 2 



Ms + Eg 



Ms — E n 



g7°g + <n°7 5 g 



{Ms - E q ) 2 - M*, I 2 L J J 

where Gp is the Fermi constant, m g is the SM quark mass, M g i is the bilepton quark mass 
and the couplings in the last two amplitudes (bilepton quark exchange) read as, g qq > = g n 
forq — d and q' = D, g qq/ = g 2 2 for q = s and q' = S and g qq > = h 33 for q = t and 
q' = T, no other quark contributes to these amplitudes. In the above equations we have used 
the non-relativistic limit, which also allows us to discard the term m(pi)7°7 5 m(p2) ~ 0. In 
this limit we can write {qj°q) ~ {qq). Besides, we can neglect the quark tri-momentum and 
write E n ~ m„. 
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As expected we will have only spin-independent (SI) contributions to the WIMP-nucleon 
cross section since our WIMP is a scalar. It is useful to define then the following average 
amplitude: 

(M) = iaJqq) , (33) 



whith a q given by, 




m q , g qq , M + m q 



M 2 H K zv PJ (M^ + m q y-M 2 q , 2 

- M ^~ m 4 (34) 

We write the matrix elements for the quarks in the nucleon by separating the contributions 
of light and heavy quarks as follows 

(qq) = ^fl? foru,d,s, 



m q 



(qq) = ^^f P T: forc,b,t. (35) 

27 m q G 

Here the subscripts p and n label the proton and the neutron, respectively. 

In order to obtain the WIMP-nucleon coupling it remains to sum over the quarks, noticing 
that all the six SM quark flavors contribute in Higgs exchange channel, while in the bilepton 
quarks exchange channels only down, strange and top quarks contribute, as can be seen from 
Eq. Q). We then get, 

fin = m p , n £ %-f%> + ^m p , n f% £ ^ , (36) 

q=u,d,s 1 q=c,b,t 1 

where = 0.020 ±0.004, = 0.026 ± 0.005, = 0.118 ±0.062, = 0.014 ± 0.003, 
f» d = 0.036 ± 0.008, = 0.118 ± 0.062 (see Ref. 0) and f*Q, which is due to the (p 
coupling to gluons through loops of heavy quarks, is obtained from, 

f P To= l - E /If> (37) 

q=u,d,s 

leading to the values, /£ « 0.84 and « 0.83. 

The WIMP-nucleus SI elastic cross section results from summing the nucleons in the 
target, yielding, at zero momentum transfer, 

with, rriN the nucleus mass, Z the atomic number and A the atomic mass. What is 
usually employed in expressing constraints on WIMP-nuclei elastic experiments is the WIMP- 
nucleon cross section, which in the SI case reads, 

2 

SI m p.n nn > 

°P,n = *0^, (39) 
where m r = M^m N j (M^ + m N ) is the reduced WIMP mass. 
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Finally, we can compute this normalized cross section for the 73 Ge (appropriate for 
CDMS spin-independent cross section) using the values of the 3-3-lRHz/ model parameters 
fixed as before. In Fig. [7] we present the current and projected data for WIMP-nucleon cross 
section from CDMS and XENON Collaborations |§, [23|, [25p, ans in Fig. § we present our 
results from 3-3-lRHz/ model. 




x x x x 

jjs ^ ifc $z $z 

+ + + + 



WIMP Mass [GeV/c ] 
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SuperCDMS (Projected) Phase B 

SuperCDMS (Projected) Phase C 



Figure 7. Spin-independent WIMP-nucleon elastic scattering cross section sensitivity (current 
and projected) from CDMS and XENON Collaborations. 



In Fig. H we have plotted only the results for V = 2 TeV and V = 3 TeV for convenience, 
where the thick regions are the allowed ones for WMAP bounds on VL C dm- However, the 
reader should have in mind that what we really have is a continuous range of WIMP mass and 
cross-section, which would fill a band between those two small regions shown in the figure. 
As we can see from Figs. [7] and ||, <\> is close to the threshold detection only for projected 
sensitivity of XENON IT experiment [E3J, but still below this threshold for the whole range 
of parameters assessed here. However, for this range of parameters, direct detection might 
be realized at improved sensitivity, lying between the projected phases B and C of Super- 
CDMS [j22"l], covering WIMP-nucleon cross-sections a little bigger than 10 -11 pb. 
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Figure 8. Spin-independent WIMP-nucleon cross sectionof the 3-3-1 RHv model as a 
function of the WIMP mass. The upper curve is for V — 3 TeV and the lower one is for 
V = 2 TeV, both taken for M q i — 1.5 TeV. The thicker regions on these lines are those in 
agreement with the bounds on Vic dm imposed by WMAP. 



In the case of neutralinos in SUSY and in 5D UED, there is also a spin-dependent 
contribution to be considered, which is dominant when A < 20, but suppressed for bigger 
values of A. For 6D UED spinless photon, similar to our case, there is only SI WIMP- 
nucleon scattering. The SI cross sections should differ basically for specific model dependent 
WIMP-quark couplings. Concerning neutralinos, typical values for this SI normalized cross 
section is a^ n « 10~ 12 - 10~ 6 pb, while for flW in 5D UED a^ n « 10~ 10 pb for WIMP 
mass about 1 TeV. In the case of spinless photon in 6D UED, a^ n ~ 10~ n — 10~ 9 pb, 
considering the region of WIMP mass for which the observed amount of CDM is generated, 
M ~ 200 — 300 GeV. In our model, only d, s and t quarks participate in WIMP-nucleon 
interaction, but this does not affect appreciably the amount of WIMP-nucleon scattering, 
yielding similar results as the UED model, namely a^ n ~ 10~ n — 10~ 10 pb for « 
600 - 1000 GeV. 

Also, collider signatures may be pursued for the forthcoming LHC at CERN |[7p, which 
will have enough center of mass energy to test several new particles at TeV scale. Collider 
signatures and indirect detection [p6|], which concerns WIMP annihilation into SM particles 
like photons and neutrinos, should be carried out in the future when more constraints on the 
parameter space of 3-3-lRHz/ model are available. Nevertheless, our aim here was to show 
that a scalar bilepton WIMP can be realized in 3-3-lRHz/ model, providing a completely 
distinct candidate for explaining CDM, reproducing its observed relic abundance and in 
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agreement with the most stringent constraints from direct detection experiments. It is true 
that we have investigated just a tiny range of the parameter space, but this limitation does not 
invalidate our conclusions and can be better explored as long as we improve our knowledge 
concerning the phenomenological consequences of this model. 



5. Conclusions 



The 3-3-lRHz/ model admits a couple of bilepton particles in its spectrum, raising the 
possibility of having a CDM candidate, since bileptons carry two units of lepton number. 
This is so because such a very specific quantum number is appropriate to forbid its interaction 
with many of the electroweak fields, since they are allowed to decay only on other bileptons, 
which have to be heavier than SM particles. Considering this scenario we obtained the particle 
mass spectrum of scalars in 3-3-lRHz/ model and, by assuming some conditions over the 
parameter space, we have shown that the lightest bilepton in the model turns out to be a 
scalar, a combination of two scalar interaction eigenstates that we called 0. For the region 
where the values of the parameters guarantee the stability of this scalar, we computed the 
abundance and obtained stringent constraints for the parameters in order to have agreement 
with WMAP results for CDM abundance. We have found that <\> can have mass ranging from 
about 600 GeV to some Few TeV, characterizing it as a heavy WIMP. It is opportune to say 
that, although we have restricted our parameter space due to lack of knowledge on several 
couplings in the model, we had no need to unnaturally adjust them to very small values as 
generally happens in several models, including Super symmetry. In fact we assumed that these 
couplings are close to one, and checked that </> is an excellent candidate to represent a WIMP 
and explain the presently observed CDM abundance. 

We also studied the possibility of observing this WIMP in direct detection experiments. 
For this we have computed the elastic scattering 0-nucleon cross section and contrasted our 
results with present and future experiments. We have seen that <p is still far from the range of 
detection for current and near future CDMS and XENON sensitivities, at least for the short 
parameter space considered in this work. However, even this limited scenario can be at reach 
for projected Phases B and C of Super-CDMS [j22l]. Besides this, it would be interesting to 
pursue the production of <p at collider experiments, mainly at LHC, and also extend our search 
including a larger region of the parameter space considering additional phenomenological 
constraints on 3-3-lRHz/ model from Collider physics and Cosmology, as well as include 
prospects for <p indirect detection too, a gap we wish to fill soon. 

Finally, we would like to stress that our proposed WIMP is not only feasible but a rea- 
sonable alternative in the sense that the Particle Physics model we are dealing with is only a 
small extension of the SM gauge group, whose scale is about to be assessed at LHC. There 
are several features that distinguishes the 3-3-lRHz/ model from other extensions, like SUSY 
and Extra dimensions models. Namely, we have not only bilepton scalars in the spectrum but 
vector bosons and quark bileptons, all of them acquiring mass at hundreds of GeV. Certainly 
their signal at detectors are worth to be studied. Besides, new phenomena are predicted in this 
model [[IT], [B], [T3J, [ZDp, including neutrinoless double beta decay, rare decays, new sources of 
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CP violation and so on. Presence of such signals would reinforce our expectation concerning 
a bilepton WIMP to explain CDM in the Universe. 
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Appendix A 

In this Appendix we show some interactions of interest for the computation of the WIMP 
abundance. 

• Interactions between WIMP and gauge bosons 

2 

^W+V^ + U.c; (A.1) 

• Interactions between WIMP and the Higgs 

- ^=(2X2V V + \6V p )H<j>*<j> - -^=(A 2 + + H.c ; (A.2) 

• Interactions between WIMP and the quarks 

- 9iad 'iLd a R<P* - h^u' ZL u aR 4> + H.c ; (A.3) 

• Interactions between Higgs and the gauge bosons 

- vVK + v p ) HW+V-» + H.c ; (A.4) 

• Interactions between Higgs bosons 
g 2 ((X 2 + X 6 ) Vri + ^-v p )HHH + H.c; (A.5) 
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